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Quantum annealing provides a promising route for the development of quantum optimization devices, but 
the usefulness of such devices will be limited in part by the range of implementable problems as dictated by 
hardware constraints. To overcome constraints imposed by restricted connectivity between qubits, a larger set 
of interactions can be approximated using minor embedding techniques whereby several physical qubits are 
used to represent a single logical qubit. However, minor embedding introduces new types of errors due to 
its approximate nature. We introduce and study quantum annealing correction schemes designed to improve 
the performance of quantum annealers in conjunction with minor embedding, thus leading to a hybrid scheme 
defined over an encoded graph. We argue that this scheme can be efficiently decoded using an energy minimiza¬ 
tion technique provided the density of errors does not exceed the per-site percolation threshold of the encoded 
graph. We test the hybrid scheme using a D-Wave Two processor on problems for which the encoded graph is a 
2-level grid and the Ising model is known to be NP-hard. The problems we consider are frustrated Ising model 
problem instances with “planted” (a priori known) solutions. Applied in conjunction with optimized energy 
penalties and decoding techniques, we find that this approach enables the quantum annealer to solve minor em¬ 
bedded instances with significantly higher success probability than it would without error correction. Our work 
demonstrates that quantum annealing correction can and should be used to improve the robustness of quantum 
annealing not only for natively embeddable problems, but also when minor embedding is used to extend the 
connectivity of physical devices. 


I. INTRODUCTION 

Quantum computers are in principle able to solve computa¬ 
tional problems that are intractable for their classical counter¬ 
parts Glia. Remarkably, rudimentary quantum annealers—a 
special type of adiabatic quantum computer O specifically 
engineered to implement quantum annealing Eii to solve 
hard optimization problems Col —are already commercially 
available ciHia, and their availability to the scientific com¬ 
munity has recently stimulated a productive debate 
The question at stake is whether current experimental quan¬ 
tum annealing can really be considered a form of quantum 
computation, even when the decoherence time of the compo¬ 
nent qubits is much smaller than the overall computational 
time. While theory suggests that the answer to this ques¬ 
tion can be positive ||25l|26]|, as it concerns quantum anneal¬ 
ing experiments the answer is still open, in particular due to 
the “black-box” nature of quantum annealing, which makes it 
difficult to unambiguously determine the nature of the physi¬ 
cal processes taking place between the initialization and read¬ 
out steps. However, convincing evidence is mounting that 
quantum effects, including entanglement IZTll and multi-qubit 
coherent tunneling 1^ . play a functional role in explain¬ 
ing the behavior of available annealers. Moreover, open sys¬ 
tem quantum dynamical models (281 |29l have been success¬ 
ful in reproducing experimental data from quantum anneal¬ 
ing experiments C6lEI]|22l|28l[30l. A closely related ques¬ 
tion is whether quantum annealers are able to provide quan¬ 
tum speedup, a theoretical possibility ED despite skepticism 


based, e.g., on the appearance of exponentially small gaps 
(3^ . While compelling evidence of a quantum speedup is 
still lacking with current physical annealers, it is also gener¬ 
ally recognized that technological improvements (in particular 
reduced control errors and shorter annealing times) in the next 
generation of quantum annealers are necessary in order to un¬ 
ambiguously address the issue ESISD. A careful choice of 
the benchmark problems is another important consideration 
(3811391, in particular in light of the fragility of spin glasses to 
small perturbations (temperature chaos) (40ll4T1l . 

Ultimately, the usefulness of quantum annealers relies on 
the promise of achieving quantum information processing on 
a large scale. This goal can only be met through the use of 
some form of quantum error correction (QEC), which is nec¬ 
essary to protect the performance of any quantum computa¬ 
tion against decoherence (42H44ll . The scalability of quan¬ 
tum computing via QEC has been established in the gate 
model in the form of the celebrated threshold theorem (e.g, 
Refs. dsma). Adiabatic quantum computation and quantum 
annealing also require some form of error correction in order 
to preserve their advantages over classical computation. This 
is true even though adiabatic dynamics is intrinsically robust 
to some forms of decoherence iiDsa. Several error sup¬ 
pression techniques suitable for adiabatic quantum computing 
(AQC) have been proposed (50H54l . However, despite valiant 
theoretical attempts (55]| much less is known about how to 
achieve fault tolerant AQC, and some negative results have 
been established (56H581 . Still, it has recently been shown ex¬ 
perimentally that it is possible to improve the performance of 
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physical quantum annealers with the help of tailored quantum 
annealing correction ISgHSTl, a method we shall employ and 
generalize in this work. 

Typically, realistic implementations of quantum annealing 
are subject to restricted (local) interactions between qubits 
specified by the physical hardware. Minor embedding (ME) 
techniques overcome this limitation by using several physi¬ 
cal qubits to represent a single logical qubit with a larger set 
of interactions |[62H65]| . A minor embedding replaces a given 
optimization problem with an approximate yet implementable 
one. As a consequence of this approximation the performance 
of a quantum annealer is typically degraded. ME also adds 
two extra steps to the computations performed with a quan¬ 
tum annealer. One is an encoding procedure that determines 
the mapping between the logical problem and its physically 
embeddable representation, which in turn requires a choice of 
energy penalties to align the physical qubits of each logical 
qubit. It turns out that performance is sensitive to the map¬ 
ping chosen and it is thus important to optimize it. The other 
is a decoding procedure that recovers the logical answer from 
the physical output. Decoding is nontrivial whenever the val¬ 
ues of the physical qubits corresponding to the same logical 
qubit disagree, and will be discussed in detail in this work. 

Quantum annealing correction (QAC) was proposed in 
Ref. 15^ , where it was demonstrated on anti-ferromagnetic 
chains and further discussed for random Ising problems in 
Ref. IbOll . These works showed that QAC can improve the 
performance of quantum annealers. Implementing several in¬ 
dependent copies of the same problem enables errors to be 
corrected using, e.g., majority vote as in classical repetition 
codes. Thus, importantly, the requirement of a purely adia¬ 
batic evolution can be relaxed since excited states can be tol¬ 
erated as long as they can be correctly decoded. Quantum 
error suppression is obtained with the use of ferromagnetic 
energy penalties connecting the various copies, that are turned 
on gradually along with the final Hamiltonian. Since these can 
also be understood as the stabilizers of the quantum repetition 
code, that detect and penalize bit-flip errors during the course 
of the evolution, agreement is enhanced within each copy by 
energetically penalizing errors that would result in disagree¬ 
ment. The combined effect of encoding into logical qubits 
and quantum error suppression can be seen as an increase of 
the effective energy scale with which the logical problem is 
physically implemented. This increases the relevant minimum 
energy gaps that protect quantum annealing from thermal and 
dynamical excitations. Moreover, it improves the response 
of quantum annealing devices to intrinsic control errors (due 
to the limited accuracy with which the physical fields can be 
tuned) by using several physical couplings to specify a single 
logical interaction EMD. 

The main goal of this work is to develop and study QAC 
techniques for improving the performance of a quantum an¬ 
nealer following the expected degradation due to the im¬ 
plementation of ME techniques. In doing so, we explore 
the efficacy and applicability of several encoding and de¬ 
coding strategies. Our tests are performed using a D-Wave 
Two (DW2) “Vesuvius” quantum annealer with 504 functional 
qubits. 


Our paper is organized as follows. In section |I^ after a brief 
review of quantum annealing, we give an overview of ME and 
QAC. We explain how ME can be viewed as a map from a 
given logical problem to the layout dictated by the physical 
implementation, while QAC reverses this direction and maps 
from the physical level to an encoding that provides protec¬ 
tion. Thus, while different in scope, the two approaches can 
be seamlessly concatenated in what we call quantum anneal¬ 
ing correction with minor embedding (QAC-ME). We also 
discuss the different role of frustration effects on ME and 
QAC. In Sec. |I^ we discuss different strategies for choos¬ 
ing the energy penalties in both ME and QAC, and with an 


eye to optimality, introduce a nonuniform strategy. In Sec. IV 


we discuss decoding strategies. These are required in order 
to assign a value to a logical or encoded qubit when the cor¬ 
responding physical qubits have different values. A majority 
vote on the physical qubits is a simple choice but that may not 
always justified. We thus propose energy minimization as a 
more effective decoding strategy and provide a general argu¬ 
ment based on percolation theory to determine whether it can 
be efficiently implemented. Section [V| introduces the “square 
code,” a QAC-ME scheme defined on a 2-level grid that can 
be implemented with the physical connectivity (“Chimera” 
graph) of the DW2 processor. The Ising problem on the 2- 
level grid is of particular interest since it was the first example 
of an NP-hard Ising model problem 1661 . The square code is 
of independent interest since it is amenable to concatenation. 
Section |Vj presents our experimental results, obtained using 
the DW2 processor. The problem instances we study are con¬ 
structed using the planted ground states method for frustrated 
loops introduced in Ref. 1341 . This method allows for a tun¬ 
able hardness parameter in the form of the effective clause 
density. The experimental results demonstrate the significant 
performance boost that can be achieved using the QAC-ME 
strategy, by exhibiting higher ground state success probabil¬ 
ities, with an improvement that is more pronounced for the 
harder problem instances. Section [Vn| reports a detailed com¬ 
parative experimental analysis of the different decoding strate¬ 
gies discussed in Sec. |lvj specialized to the square code. In 
Sec. VIII| we present the results of a simulated quantum an¬ 
nealing study that numerically investigates the role of tem¬ 
perature and control noise on the performance of the ME and 
QAC-ME strategies, as these are parameters that we cannot 
control experimentally. We conclude and discuss possible de¬ 
velopments of our work in Sec.|I^ Additional supporting in¬ 
formation is presented in the Appendix. 


II. QUANTUM ANNEALING, MINOR EMBEDDING, AND 
QUANTUM ANNEALING CORRECTION 

A. Quantum annealing 

The term quantum annealing was first coined (as far as we 
know) in 1994 by Einnila et al. (H, though the underlying 
idea was proposed at least as early as 1984 Ezl as a numeri¬ 
cal algorithm inspired by simulated annealing 1^ . where the 
role of thermal fluctuations was conjectured to be advanta- 
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(a) USC-ISI DW2 Chimera graph (b) Annealing schedule 

FIG. 1. (a) Physical connectivity of the “Chimera” graph of the 512-qubit DW2 processor. White circles correspond to the 8 unusable qubits in 
the USC-ISI DW2 processor, (b) Annealing schedule of the DW2 processor. The dashed horizontal line is the physical temperature of 17mK 
at which the device is kept. 


geously replaced by quantum tunneling. This was soon fol¬ 
lowed by numerous other studies based on a similar idea (see 
Refs. |[5l [69ti7^ for some of the earliest works and Ref. (731 
for a recent review). Numerical and experimental results have 
suggested that quantum tunneling can be more effective than 
thermal fluctuations for reaching the ground state iH). This 
conclusion has recently been revisited CD, and remains the 
subject of active scrutiny (e.g., Ref. (751 ). 

Quantum annealing for spin systems is typically achieved 
by the following time-dependent, transverse-fleld Ising 
Hamiltonian: 

H(t) = A{t)Hx + ^ e [0, ty] . (1) 

The term Hx = (with erf being the Pauli matrix act¬ 

ing on qubit i) is a transverse fleld whose amplitude controls 
the tunneling rate. The solution to an optimization problem of 
interest is encoded in the ground state of the problem Hamil¬ 
tonian Hp, given by 

Hp=J2 hi<Jt + Y. ’ (2) 

where the sums run over the vertices V and edges £ of the 
weighted, undirected graph G = {Vh,£j), where the local 
flelds hi and couplings Jij are the weights. In a closed system, 
the adiabatic theorem ensures that if the system is initialized in 
the ground state of 77(0) = A{0)Hx (assuming B{0) = 0), 
a sufficiently slow (adiabatic) evolution results at the end of 
the computation in the ground state of the final Hamiltonian 
H{tf) = B{tf)Hp (assuming 74(ty) = 0). In an open system 


annealing at a non-zero temperature one expects the system to 
evolve to a mixed final state, e.g., the Gibbs state of 77(ty), 
if equilibration is possible throughout the annealing process 
(see, e.g.. Ref. fISl ). 

The D-Wave devices use an array of superconducting 
niobium-based flux qubits in order to physically realize 
quantum annealing (TTHT3]| . Each vertex of the D-Wave 
“Chimera” graph is a physical flux qubit and each edge is 
a physical inductive coupling between two qubits. Fig. |l(a)| 
shows a pictorial representation of such a graph. In the DW2 
processor used in this work the temperature is set to 17mK 
and the annealing sche dule i s given by the functions A{t) 
and B(t) shown in Fig. 1(b) In the following, we specify 
the values of the couplings hi and Jij in dimensionless units, 
with the maximum strength of the allowed couplings given by 
|^max| ^ 2and = 1. 


B. Mapping to and from the logical problem: ME and QAC 

The problem Hamiltonian defined in Eq. © may or may 
not be specified in terms that directly correspond to physical 
variables representing the actual degrees of freedom of a phys¬ 
ical device. When such a specification is possible we have a 
direct (D) embedding. An example is the class of random 
Ising model problems studied in Refs. EllSl, which simply 
require a specification of the hi and Jij on the Chimera graph 
of the D-Wave device. 

However, often one is interested in solving problems that 
do not have such a direct embedding. E.g., consider a prob- 
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FIG. 2. ME: Schematic of two logical groups i, j (large circles) rep¬ 
resenting two logical qubits a-|. Physical qubits are denoted by 
blue circles; physical problem couplings are denoted by black lines, 
penalty couplings by red lines. In the example shown the penalty 
coupling and the physical problem couplings Jj are 

active, among others. 


lem defined on a complete graph, which requires a minor em¬ 
bedding on a graph that is not complete. In addition one might 
want an encoding of the Hamiltonian in order to perform error 
correction. To deal with this more general scenario we shall 
proceed to define two maps: one for ME, another for QAC. 
We shall also consider the concatenated map, of ME followed 
by QAC, which we refer to as QAC-ME. 

7. ME: from the logical Hamiltonian to the physical Hamiltonian 

When the problem Hamiltonian is defined in terms of vari¬ 
ables that do not have a direct representation and requires a 
minor embedding we shall refer to it as a “logical Hamilto¬ 
nian” and write 

= E Wl- (3) 

(I3)e£j 

Quantities decorated with a tilde represent the abstract logi¬ 
cal local fields, couplings, Pauli operators, etc. The logical 
Hamiltonian is defined on the vertices and edges of the logi¬ 
cal graph G = (V^, Ej). Given the logical Hamiltonian which 
specifies the optimization problem of interest in terms of ab¬ 
stract variables, the task is to map it to a physical problem 
Hamiltonian that has a direct embedding. This is the minor 
embedding map 

Hp Hp . (4) 

This map exists when (5 is a minor of G, i.e., it has a minor 
embedding in gQ 


^ A minor graph is obtained from G by collapsing groups of connected 
(physical) vertices of G into single (logical) vertices of G and possibly 
removing edges from gES). Note that TWme is a multivalued map (one- 
to-many). 



FIG. 3. Schematic example of the map AIme- Two logical qubits 
i and j (red dots) with four logical couplings each (black lines) are 
expanded into two physical qubits each (circled blue dots) with three 
physical couplings each. The red lines connecting the physical qubits 
are the penalty couplings and • In this example the logical 
coupling J~fj is replaced by a single physical coupling special 

case of Eq. 


Each of the TV = |V^| vertices of G is occupied by a clas¬ 
sical binary (±1) variable, which we call a logical qubit|^ 
The map A^me replaces these by TV = |V/i| physical qubits. 
A minor embedding defines a partitioning of all TV physical 
qubits into TV disjoint logical groups {i}, each comprising 
N~- physical qubits, whose union covers all physical qubits, 

i.e, TV = Note that we distinguish between a log¬ 

ical group and a logical qubit. The latter was defined above; 
the former is a set of physical qubits representing the logical 
qubit. What identifies a logical group is a set of ferromagnetic 
“penalty” couplings which tie the physical qubits representing 
a given logical qubit together, as illustrated in Eig. Penalty 
couplings are always “intra”, i.e., couple physical qubits of a 
given logical group i, and can thus always be written as 

Jpen ^ 

The physical qubits comprising two logical groups are cou¬ 
pled in a manner that depends on the problem Hamiltonian. 
The logical local fields and couplings are defined via 

N-i N-i ^3 

k = k = T.T,\i, « s ) 

k=l k=l 1=1 

where some of the may vanish, as illustrated in Eig. l2j 
Note that Eq. ([^ does not uniquely specify the physical locm 


^ Technically this is a misnomer since the variable is classical. However, we 
have in mind the entire time-dependent Hamiltonian of quantum annealing, 
and moreover, we can always measure the logical qubit in the logical 
basis. 
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fields and couplings (which are always “inter”), allowing for 
a large arbitrariness in fixing the mapping A^me between the 
logical and physical problems. This arbitrariness is in fact 
useful, since the efficacy of an ME implementation can be 
very sensitive to the specific choice of the physical couplings. 

Fig. shows schematically how a logical-graph vertex of 
degree four can be minor embedded into two physical-graph 
vertices of degree three, and how two logical qubits are corre¬ 
spondingly replaced by four physical qubits. 

Ideally, a minor embedding should give a faithful repre¬ 
sentation of the original logical problem, in the sense that it 
should encode the solution to the logical problem in its ground 
state. Intuitively, this can be achieved by implementing large 
penalty couplings, though it is clear that a balance must be 
struck so as not to overwhelm the problem couplings. Con¬ 
sequently the optimal choice for the energy penalties is not 
obvious. Very large penalties make logical groups behave 
as independent ferromagnetic clusters. These penalties can 
overwhelm the relatively weaker interactions between differ¬ 
ent clusters that represent the logical problem. When the log¬ 
ical groups are small, theoretical and experimental analyses 
suggest that the optimal value of the energy penalties is of 
the same order of magnitude as that of the logical couplings 


i j 


FIG. 4. QAC: Schematic of two encoded groups i,j (large circles) 
representing two encoded qubits af, . Physical qubits are denoted 

by blue circles; physical problem couplings are denoted by black 
lines, penalty couplings by red lines In the QAC case the problem 
couplings are an injection (one-to-one but not necessarily onto) be¬ 
tween the physical qubits of each encoded group. As depicted, the 
penalty couplings may involve all or a subset of all physical qubits 
of each encoded group, and may vary from group to group. 


tonian” and write 








(7) 


Hp = j2 hn + E • (8) 


We discuss how this freedom can be exploited in order to im¬ 
prove the performance of ME implementations in Section [nl| 


An important consequence of the finiteness of the energy 
penalties within each logical group is that, informally, logi¬ 
cal groups can “break”, i.e., the values of the physical qubits 
(measured in the eigenbasis at tj) representing a logical 
qubit can disagree. This may result in an unfaithful represen¬ 
tation of the logical qubit by its logical group. We call logi¬ 
cal groups with internal disagreement “broken qubits” (BQs). 
More formally, the logical qubit belongs to the code subspace 
defined by the two “all-agree” states of the logical group, and 
any “disagreement” is a non-code state (more on this in our 
discussion of QAC below). Because of Eq. 0 BQs are ubiq¬ 
uitous in every ME implementation. Classical post-processing 
(decoding) strategies that assign appropriate logical values to 
BQs are thus a necessary aspect of minor embedded quantum 


annealing. We discuss several such strategies in Section IV 


2. QAC: from the physical Hamiltonian to the encoded 
Hamiltonian 

When the problem Hamiltonian Hp in Eq. ^ needs to be 
encoded in order to provide protection against decoherence 
and noise we shall refer to the result as the “encoded Hamil- 


^ When long chains are used as logical groups, the optimal value for the 
energy penalties seems to scale linearly with the chain length Ell. 


Quantities decorated with a bar represent the encoded opera¬ 
tors, acting on the N encoded qubits. The encoded Hamilto¬ 
nian is defined on the vertices and edges of the encoded graph 
G = (V/i, £j). Note that unlike the ME case, the same local 
fields and couplings are used in the encoded Hamiltonian as 
in the physical Hamiltonian (i.e., h and J do not have a bar), 
except that they are defined on the encoded graph. In particu¬ 
lar, h-^ = hi and J-fj = Jij • Given the physical Hamiltonian, 
the task is to map it to an encoded Hamiltonian that also has 
a direct embedding. This is the quantum annealing correction 
map 


Hp Hp . (9) 

This map exists when G is a subgraph of G. I.e., QAC re¬ 
quires several copies of the logical graph G to be directly em¬ 
beddable into the physical graph G. 

Each of the = \Vh\ vertices of G is occupied by a clas¬ 
sical binary (±1) variable. The map Mqac replaces these by 
N = \Vh \ encoded qubits. Similarly to the ME case, we iden¬ 
tify encoded groups as sets of physical qubits tied together by 
ferromagnetic penalty couplings, as depicted in Fig.|^ Unlike 
the ME case, we impose that the problem couplings between 
encoded groups are injective (one-to-one) and that all encoded 
groups are of identical size: = AT Vi G {!,...,7V}. The 

penalty couplings may be non-uniform in order to allow for 
their optimization, e.g., as in Eq. Q. 

In the QAC case encoding acts to provide protection. To see 
how this works, we briefiy digress to frame our discussion in a 
more general setting. Let us recall that one can define a quan¬ 
tum error correction code via a set of stabilizers, {Sg}, (or 








6 


parity checks in classical error correction language) 14411781 . 
This set of commuting operators on n qubits defines a col¬ 
lection of syndrome subspaces, each encoding a fixed num¬ 
ber of encoded qubits k with corresponding encoded opera¬ 
tors (we temporarily drop the bar from the qubit 

subscripts; we will sometimes return to this notation if there is 
no danger of confusion). These subspaces are defined accord¬ 
ing to their eigenvalues with respect to the stabilizer operators. 
Since our focus is on encoding the final Hamiltonian we 
only consider bit flip errors and thus restrict our attention to 
quantum repetition codes. 

An [[n^k^d]] quantum repetition code (encoding n physical 
qubits into k encoded qubits and having distanced = n can 
be defined by the stabilizers crfcr| where {i, j} can be cho¬ 
sen as all the nearest neighbor pairs. Such a code can correct 
t = [d/2\ errors, and can be reliably decoded using majority 
voting if the number of errors does not exceed t. The decoding 
is undecided in the case of a tie, i.e., for n/2 errors (n even) 
and fails if the number of errors exceeds t. 

For simplicity we consider the case of a single encoded 
qubit {k = 1) per encoded group of n, and consider multi¬ 
ple encoded groups {t\. Such a code has the basis states 

|0) = |0)®” ; |I) = |1)®" . (10) 

The encoded operators are = (cr“)®" and, e.g., = erf. 

The QAC protection idea is to exploit the various realiza¬ 
tions of the encoded operators in order to provide an energy 
boost, and to include stabilizer terms in such a way that errors 
causing transitions outside of the code space are penalized. In 
particular, writing the unencoded final Hamiltonian [Eq. 0] 
in the form 

hial + y] Ju,(jI(jI, , (11) 

I e,e' 

the encoded final Hamiltonian becomes 

H = ^hiY^[(7^]r + ^^Ju’[(Te<^e]r , (12) 

i r r 

where [A]r is the rth realization of the operator A. These 
multiple realizations exist since the logical Pauli operators are 
equivalent under multiplication by the code stabilizers, and for 
our repetition code any tensor product of an odd number of 
operators is an equivalent realization of . This multiplicity 
boosts the energy scale of H relative to H, and it increases 
the gap against thermal excitations. This is complemented by 
adding to ^ a penalty term that is a sum over the available 
stabilizers Sa for each encoded qubit i\ 

J^pen = -Jpeny][5.]«- (B) 

i,s 


^ A code can detect any Pauli errors error E having weight < d — 1, i.e., 
(V^al E I'lpi)) = eSab for any pair of codewords \'ipa) , I'lpb) and constant 
c. Equivalently, the distance is the minimum weight of a Pauli operator E 
such that {'ipa\ E I'lpb) ^ c5ab- 


Each bit-flip error anticommutes with at least one term in 
i7pen and hence incurs an energy cost of at least 2 Jpen- Such 
errors are thus thermally suppressed (sa [531 El ESI [571 [59]- 

ISB. 

Returning now to our specific setting and previous notation, 
note that whereas in the ME case the were abstract binary 
variables, in the QAC case we choose the encoded Pauli oper¬ 
ators as the total Pauli operators of each encoded qubit, i.e., 

K' 

= (14) 

k=l 

where we allow for the possibility that not all K physical 
qubits in each encoded group participate in defining the en¬ 
coded Pauli operators (i.e., K' < K). In this manner the 
encoded problem Hamiltonian is expressed in terms of the en¬ 
coded operators of the quantum repetition code. This means 
that 1 -^ J-A-Ad^d-. = J-a-a a-. , i.e., each cou- 

^3 I 3 ^3 i 3 ^3 ^k=l Ik 3k' 

pling Jij is boosted by a factor of K', since on the encoded 
graph each J-q is represented by K' coupled pairs of physical 
qubits. We stress again that this boost is a useful feature since 
it increases the energy scale of the encoded problem Hamilto¬ 
nian, a feature that works to suppress both thermal excitations 
and random control errors EaMD- 

Broken encoded qubits appear in the QAC case just as BQs 
appear in the ME case, manifesting as encoded groups whose 
physical qubits yield different results when measured at t/, 
resulting in a measured eigenvalue of df that is smaller than 
K' in absolute value. Our QAC code can be viewed as a 
1, K']] repetition code, i.e., a code with distance K' en¬ 
coding one qubit into K'. The penalty terms J-a. j.a? im- 
plicit in Hp can be understood as elements of the stabilizers 
of the repetition code. Since these stabilizers detect bit-flip 
errors, as mentioned above these errors are energetically sup¬ 
pressed. 

The transverse field is invariant under the transformation to 
the encoded operator as in Eq. (T4\ \ 

N N K N 

Hx = = EE< = E^f = Hx ■ (15) 

*=1 i=l k=l i=l 

This observation is motivated by the fact that the D-Wave pro¬ 
cessors do not allow control over the transverse field. How¬ 
ever, we note that for the purposes of error correction this rep¬ 
resents a significant limitation, since (as mentioned above) 
the encoded operator of the repetition code is in fact 
d^ = (8)/cCr^, and the inability to implement this (many-body) 
operator compromises the efficacy of QAC lEHEOl as well as 
other strategies 1501 [561 . 


C. ME vs QAC: frustration effects 

QAC is known to be effective in improving the performance 
of a quantum annealer when compared to a direct (D) embed¬ 
ding of optimization problems |59l[60l. On the other hand. 
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FIG. 5. Contrasting ME and QAC via frustration effects. ME: The 
top configuration displays an example of a frustrated logical qubit 
(red) with ferromagnetic logical interactions (black). Erustration is a 
consequence of the opposite orientations of the outer logical qubits 
(represented as numbered arrows). Erustration of the logical interac¬ 
tions can be removed by the appearance of a broken qubit (lower 
left) after minor embedding. This configuration can be energeti¬ 
cally favored if the penalty coupling (red line) is sufficiently weak. 
Thus, in the ME case, frustration may favor the appearance of bro¬ 
ken qubits. QAC: Now the top configuration represents a frustrated 
physical qubit (red) with ferromagnetic couplings (black). In this 
example the QAC map replaces this configuration with two copies 
(bottom right), coupled with a ferromagnetic penalty (red). After the 
QAC map it is energetically favorable for the penalty-coupled qubits 
(blue) to be aligned for arbitrarily small values of the penalty. Thus, 
in the QAC case, frustration does not favor the appearance of broken 
qubits. Instead, configurations with broken qubits are excited states 
of the encoded problem for arbitrarily small values of the penalty 
terms. 


ME typically reduces performance relative to a direct imple¬ 
mentation of a given optimization problem. This opposite re¬ 
sponse to ME and QAC has an intuitive explanation in how 
ME and QAC behave under frustration, which can be use¬ 
ful once the transverse field is off and the evolution is ef¬ 
fectively classical [i .e., w hen A{t) kT < B{t), e.g., for 
t/tf ^ 0.8 in Eig. l(b)|. The problem Hamiltonian corre¬ 
sponding to a hard optimization problem typically represents 
a strongly frustrated spin system, i.e., wherein the ground state 
configuration cannot satisfy all logical couplings. An exam¬ 
ple of what happens when a frustrated logical qubit is minor 
embedded or QAC embedded is given in Eig. In the ME 
case (where we assume that the logical problem cannot be 
directly embedded), if the penalty term is sufficiently weak, 
frustration results in a broken qubit, thus introducing new log¬ 
ical errors (to-be-decoded broken qubits). In the QAC case 
(where we assume that the problem can be directly embed¬ 
ded), in contrast, each copy of the problem bears the same 
frustration independently, and it is energetically favorable to 
satisfy the penalty coupling, no matter how weak it is. There¬ 
fore in QAC frustration is not the reason for the appearance 
of logical errors (whose source can instead be, e.g., thermal). 



EIG. 6. Schematic comparing the performance of different embed¬ 
ding schemes. The time to solution (TTS) axis is not to be taken lit¬ 
erally and represents wishful thinking. The goal of QAC is to reduce 
TTS (or increase the success probability) with respect to a D embed¬ 
ding, while QAC-ME is meant to improve over ME. Due to various 
tradeoffs it is not generally clear whether QAC-ME improves upon 
D. Of course, a D embedding may not be available, in which case 
QAC-ME is the only option. 


Rather, thanks to their connection via energy penalties, the 
different copies of the problem help each other to relax more 
efficiently to the logical ground state. 


D. Quantum annealing correction on minor embeddings 
(QAC-ME) 


It is natural to consider a concatenation of the ME and QAC 
maps: 


Hp Hp Hp . (16) 

This is the QAC-ME map. It replaces a logical problem 
Hamiltonian by an encoded problem Hamiltonian, after a mi¬ 
nor embedding of the former. 

The ME and QAC maps both expand the number of qubits. 
In the ME case a logical qubit is expanded into several physi¬ 
cal qubits in order to enable an embedding of the logical prob¬ 
lem Hamiltonian into a physical problem Hamiltonian that 
obeys the constraints of a given hardware graph. In the QAC 
case a physical qubit is expanded into several physical qubits 
in order to generate protection against decoherence and noise. 

Eigurej^ summarizes the various embeddings we have dis¬ 
cussed so far. A problem formulated on a logical graph G 
can be directly embedded (D) on a physical device with given 
hardware graph G if (5 is a subgraph of G. QAC aims to im¬ 
prove the success probability of solving this problem as com¬ 
pared to a D implementation. On the other hand, for problems 
formulated on a logical graph that can only be minor embed¬ 
ded, the goal of QAC-ME is to improve the success with re¬ 
spect to ME. It is typically not possible to directly compare 
D/QAC and ME/QAC-ME, simply because direct implemen¬ 
tations of a given optimization problem are generically un¬ 
available. It is however of great interest to consider cases 
where both D and ME/QAC-ME implementations are possi- 
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ble. This gives us an estimate of what performance decay is 
due to ME , and how much of it can be recovered by quantum 
annealing correction schemes such as QAC-ME. We perform 
such comparisons below. 


III. CHOOSING THE PENALTIES 

The penalty strengths should be optimized in order to max¬ 
imize the performance of the embedding procedures. It is 
impractical to perform this optimization on an instance-by- 
instance basis, so it is important to find a general rule for de¬ 
termining the optimal choice of penalties. This is an open 
problem that we address below. 


A. Homogeneous penalties 

Recent work has shown both experimentally and theo¬ 
retically that even the simplest choice—a constant—for the 
penalties, 

'^pen — T 5 

allows for successful QAC 15^1^ or ME fTTll . The optimal 
value of 7 is found to be roughly proportional to the value 
of the physical couplings that represent logical connections 
and typically grows with the size of the logical (or encoded) 
groups. 


B. Nonuniform penalties 


IV. DECODING STRATEGIES 

A decoding strategy is a rule that assigns a value to each 
logical or encoded qubit given the values of the physical 
qubits retrieved from the annealer in the read-out process. As 
already mentioned in the general discussion about ME and 
QAC, in the decoding process we first distinguish between 
“broken” and “unbroken” logical or encoded qubits. The un¬ 
broken case corresponds to a logical or encoded group of 
physical qubits that all have the same measured final values. 
In this case, decoding is trivial and consists in assigning to 
the logical or encoded qubit the common value of the corre¬ 
sponding physical qubits. The broken case, on the other hand, 
corresponds to a logical or encoded group where there is at 
least one physical qubits whose value differs from the others. 
Decoding strictly applies to BQs, and it can be done in dif¬ 
ferent ways. We will consider two general approaches (and a 
combination therefore). The first one is “local”, in the sense 
that the decoding assignment is performed on each logical or 
encoded group independently. It is the usual strategy for de¬ 
coding classical repetition codes and relies on majority vote 
and coin tossing. The other can be considered a “global” ap¬ 
proach, in which the value of a logical qubit is decoded by 
a global energy-minimization procedure that simultaneously 
involves all BQs. 


A. Coin tossing (CT) 

The simplest decoding strategy is random decoding (“coin 
tossing”), wherein we assign a random ±1 value to each logi¬ 
cal or encoded qubit. This strategy merely serves as a baseline 
against which we compare all other strategies. 


The optimal choice of penalty strengths achieves a delicate 
balance of tying the cluster of physical qubits forming a logi¬ 
cal group or encoded group together while not overwhelming 
the problem couplings. We can expect that a better balance 
can be achieved with a choice of penalties that depends on the 
local structure of the logical problem. One possibility is to 
choose the energy penalties as follows: 

ME: Vfc,; 4 = mean^lJijI (18a) 

QAC. VA:, / = ^pen,2 T I ,jf I 5 (1^^) 


These expressions assign a value to the energy penalties as¬ 
sociated with logical qubit i (encoded qubit i) that is propor¬ 
tional to the mean of the absolute strength of the logical cou¬ 
plings connected to that logical qubit (encoded qubit). The 
overall constant factor 7 is then optimized as in the uniform 
case for the best results. Typically, practical problems trans¬ 
late into Ising Hamiltonians with a certain amount of struc¬ 
ture, both in the connectivity and in the value of the logical 
couplings. We thus expect this nonuniform choice to be gen¬ 
erally more effective. In section VI we present a study of the 
efficacy of the uniform and nonuniform encoding strategies 
described above using the DW2 device. 


B. Majority vote (MV) 

In the majority vote strategy the value of each encoded 
qubit is chosen as the most recurrent value within an encoded 
group: 


Si = sign(y] s-J , (19) 

k 

where s and s stand for the values of the encoded and physical 
qubits in the basis ait = tf, respectively|^MV decoding is 
most successful under the assumption of independence of the 
error probability Perr for each physical bit. A repetition code 
of length n can be reliably decoded with up to d = \n/ 2 J er¬ 
rors, an event that occurs with probability which is expo¬ 
nentially decreasing in the code distance. An important caveat 
of applying MV to ME and QAC is that physical excitations 
above the ground state occurring during the evolution cannot 


^ From hereon, if there is no danger of confusion, we use the QAC terminol¬ 
ogy and notation to denote both QAC and ME. 
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FIG. 7. Failure modes for majority vote in ME. A logical qubit is 
represented here by a physical chain. MV decoding on the bottom 
chain gives a spin-down result, on the top it gives spin-up. The chain 
can break iso-energetically in different places and decoding by ma¬ 
jority vote, in this case, is likely to give no better results than coin 
tossing. 


be considered as generating independent errors on the physi¬ 
cal qubits. MV itself thus does not ensure exponentially im¬ 
proving resiliency of the ME implementation as a function of 
the size of the encoded group. 

In QAC schemes, where the encoded qubits are imple¬ 
mented as several equivalent physical copies, MV will be suc¬ 
cessful if errors in the different copies can be considered as in¬ 
dependent. A quantum contribution in the QAC scheme arises 
from the energy penalties, which provide further suppression 
of errors. MV has indeed proven to be an effective approach 
to decoding QAC IS^IbOll . 

However, MV is problematic for most ME schemes. In this 
case the absence of BQs is crucial for the correct represen¬ 
tation of the logical problem as a minor embedded problem, 
since a logical qubit can break in several different ways, giv¬ 
ing rise to different MV-decoded logical values. These differ¬ 
ent “breaking modes” may have similar energy cost, making 
them similarly probable. A typical example of this is shown in 
Eig.[7j where a logical qubit is represented by a physical chain. 
The chain can break in several places, all of them costing the 
same energy, giving different results through MV. We should 
not expect, in this simple case, MV to be more effective than 
decoding by coin tossing. 

Another issue that arises for MV is how to decode a logi¬ 
cal qubit comprising an even number n of physical qubits. A 
tie (n/2 errors) cannot be decoded by MV, or by any “local” 
decision rule other than coin tossing. 


C. Energy minimization (EM) 

Since we are interested in optimization problems, a natural 
alternative to MV is a “global” decoding strategy wherein all 
BQs are decoded together in such a way as to minimize the en¬ 
ergy of the encoded qubits. Decoding a particular state with 
this energy minimization strategy corresponds to solving an 
optimization problem on the subgraph of the encoded graph 
induced by the set of BQs, i.e., the qubits that require decod¬ 
ing. This problem is equivalent to minimizing the following 


“decoding” Hamiltonian 

^ ^ Jj-jS-j I ^f+ ^ • (20) 

zGBQ \ jVBQ / ijGBQ 

This energy minimization strategy is rather general and can be 
used to decode ME, QAC or QAC-MEC schemes. Moreover, 
it is applicable to cases where logical qubits comprise an even 
number of physical qubits, since it can also resolve ties. 

A disadvantage of EM is that, since it requires solving a 
global optimization problem, it can be computationally ex¬ 
pensive. If the number of BQs is small (in a sense we quan¬ 
tify below), this optimization can be done exactly; otherwise, 
heuristic methods are necessary. We discuss EM decoding via 
simulated annealing, along with experimental results, in Ap¬ 
pendix In the next section we argue that for a sufficiently 
small number of BQs, EM decoding can be implemented effi¬ 
ciently. 


D. Efficient optimal decoding and the percolation threshold 

We now give a general criterion to determine whether en¬ 
ergy minimization decoding can be done efficiently. We use 
well-known results of percolation theory (Ill and rely on the 
simplifying assumption that an encoded qubit is broken with a 
qubit-independent probability Pbq- This probability depends 
on the strength of the energy penalties used for the encoding. 
Eventually, the relevant value of pbq is obtained when it is 
evaluated at the optimal penalty value. 

Let us start from two extreme cases. Eirst, if Pbq ^ 1 er¬ 
rors are sparse, and it is unlikely that two BQs are coupled. In 
this case the decoding effort grows only linearly with the size 
of the problem because each BQ can be optimized indepen¬ 
dently. In such a regime, the EM strategy can be implemented 
efficiently, since the decoding time scales linearly in the num¬ 
ber of encoded qubits 


Tj^{N)r^N. ( 21 ) 

Second, if pbq is close to 1, most BQs are likely to be cou¬ 
pled. In this regime the EM decoding strategy requires one to 
solve an optimization problem that is as hard as the original 
problem: 


TD(iV)-Tp(7V), (22) 

where Tp is the time to solve the original problem. In this 
regime, EM decoding cannot be implemented efficiently. 

In general, percolation theory can be used to establish a 
“decoding threshold” pd, such that EM decoding can be im¬ 
plemented efficiently (in the large N limit) if and only if 
Pbq < Pd [this assumes that pbq is roughly constant across 
the graph, an assumption we test and confirm experimentally 
below (see Eig.[^]. It turns out that pd is equal to the “per- 
site” percolation threshold pth of the encoded graph. To see 
this, we note that the decoding hamiltonian is defined on 
a “per-site” percolation subgraph, (i.e., the subgraph induced 









10 


by vertices selected with some probability p). A fundamental 
result of percolation theory is that (for TV 1 ) there exists 
a threshold probability pth, determined by the local structure 
of the graph, above which an infinite (i.e., spanning) cluster is 
present with probability 1. This result implies that EM decod¬ 
ing cannot be done efficiently if pbq > Pth, since it requires 
minimizing for a system of BQs that is as large as the 
original problem. 

On the other hand, below the threshold pbq < Pth, the 
probability to have an infinitely large cluster vanishes. More 
precisely it can be shown that the probability that a given BQ 
is connected to at least f other BQs decays exponentially with 
the size f of the connected domain ISOl . Assuming that the 
N encoded qubits grow their domains independently, we can 
thus express the total probability to have no domains of more 
than f connected encoded qubits as 


P{r) = {1-ae ^1-aNe , {pbq < Pth) 

(23) 

for some constants a and 7 and for a sufficiently small expo¬ 
nential factor. Let us now fix the largest size of a decodable 
domain to be f. If the encoded state contains a domain of 
larger size, we deem it “undecodable,” and we reject it. An 
undecodable event will happen with probability 1 — P(f). Let 
us again denote by Tp(f) the time needed to optimize over 
a domain of size f (assuming that the decoding optimization 
is at most as hard as the original optimization problem). The 
total time required for decoding Td (TV) is proportional to the 
typical number y of independent domains of size f in the sys¬ 
tem, multiplied by the optimization time Tp (f) of a region of 
size f. Also, we note that the higher the probability of hav¬ 
ing no domains of more than r connected encoded qubits, the 
shorter is the total decoding time. Thus: 


Td(TV) - 


N Tp(f) 
f P(f) 


(24) 


Equation shows that if we take r ^ log(TV), P(r) can be 
kept constant and close to 1 as TV ^ cx). Thus we have 


TV 

?b(A^)~ j^^Tp(log(7V)), (25) 

which is the main result of this section. Below the percola¬ 
tion threshold, as long as we do not attempt to decode do¬ 
mains whose size scales faster than log (TV), the time required 
for optimal decoding through energy minimization scales sub- 
exponentially as compared to the time to find a solution to the 
original problem. In the large size limit, therefore, the decod¬ 
ing effort is negligible. 

We note that, in principle, the constant prefactor needed 
to have a small P(f) can be very large and it diverges when 
we approach the percolation threshold. However, because the 
percolation transition is usually very sharp, we can expect the 
prefactor to be small if we are well below the threshold. 

We revisit the percolation threshold experimentally below, 
and show (see Lig.[^ that for the problem instances we con¬ 
sidered in this work we are well below the threshold, so that 


efficient decoding is possible. 


V. QAC-ME OF AN NP-HARD PROBLEM 

In this section we describe an explicit QAC-ME scheme 
that combines a well-known NP-hard problem with a new 
QAC code. 


A. The two-level-grid (2LG) graph and its embedding 


Consider a logical graph G consisting of two connected 
square lattices, i.e., a two-level-grid (2LG) as depicted in 


Eig. | 8 (a)| The Ising spin glass problem on this graph is of 
interest, among other reasons, since it was perhaps the first 
shown to be NP-hard, for couplings in { — 1, 0,1} and no local 
fields (Barahona’s problem P3 ll 66 l ). This problem does not 
have a direct embedding in the Chimera graph, but it can be 
minor-embedded, as illustrated in Pig. | 8 (b) 


The ME employs logical groups with pairs of physical 
qubits connected by a ferromagnetic coupling. Logical cou¬ 
plings are represented by one or two physical couplings per 
physical qubit. Note that fewer than half of the physical cou¬ 
plings are used in this ME. The remaining physical couplings 
can be used to implement a QAC-ME scheme as shown in 
Eig. | 8 (c) using a new QAC “square” code we describe below. 


B. Direct embedding of the 2LG 


Since we anticipate that minor embedding will perform 
poorly when compared to a direct embedding, we are also in¬ 
terested in studying a direct embedding of the 2LG. An inter¬ 
esting set of such problems includes all instances that can be 
defined on the subgraph of the encoded graph of the square 
code. An example is shown in Eig. |9(a)| which depicts the 
largest 2LG subgraph directly embeddable into the USC-ISI 
DW2 Chimera graph. An example of a direct embedding of 
this subgraph is shown in Eig. |9(b)| 


C. The square code and concatenation 


We now describe the square code in detail. Using the nota¬ 
tion [[n, /c, d]]p to denote a distance d code that uses n physical 
data qubits and p additional penalty qubits to encode k logical 
qubits, this code can be understood as a [[4,1,4 ]]q codeH 
The square code has a direct embedding into the Chimera 


graph, as shown in Eig. 10(a)| This yields a 2LG encoded 
graph, as shown in Eig. |10(b)T that can be extended as shown 
in Eig. | 8 (a)[ with extent bounded only by the size of the phys¬ 
ical Chimera graph. 


^ The code introduced and used in Refs. is a [[3,1, 3]] code in this 

notation. 
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(a) 2LG logical graph connectivity 



FIG. 8. From a logical problem to a QAC-ME embedding, (a) Barahona’s “two-level-grid” (2LG) graph, for which the Ising spin glass problem 
with couplings in {—1, 0,1} is NP-hard 1661 . Disconnected vertices indicate spins with couplings set to zero, as well as unused logical qubits 
in the logical DW2 Chimera graph, (b) Minor embedding of the 2LG graph into the physical DW2 Chimera graph. White circles correspond 
to unused or unusable qubits, (c) QAC-ME embedding of the 2LG problem using the “square” code. In (b) and (c) penalties are represented 
by red (thick) couplings between groups of two (ME) and four (QAC-ME) physical qubits. Black (thin) links implement logical couplings. 
See also Eig.[^for a more detailed description of the QAC-ME encoding. 
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EIG. 9. (a) The largest subgraph of the 2LG (encoded graph of the square code) that can be directly embedded into the USC-ISIDW2 Chimera 
graph. To see that this is a subgraph, note that some of the edges depicted in Eig. |8(a)| are missing, e.g., the edge between vertices 1 and 2. (b) 
Direct embedding of (a) into the Chimera graph. Each physical qubit in (b) corresponds to a logical qubit in (a). 


As indicated by the pairs of orange and green lines connect¬ 
ing pairs of encoded qubits in Fig. 1 10(a)] each encoded qubit 
can couple to its neighbor via a 0 coupling in two dis¬ 
tinct ways, resulting in a boost by a factor of 2 at most. On 
the other hand, there are 8 ways of realizing a blue logical 


0 coupling, but, since we wish to maintain the relative 
strength between qubit-qubit coupling terms, we use only 2 of 
these realizations, or alternatively we can use all 8 but with 
1/4 the value of J per physical coupling. Similarly, there are 
4 ways of realizing a logical operator by addressing each 
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(a) Arrangement of physical qubits 
in the Chimera graph 


(b) Section of the 
corresponding encoded 
graph 


FIG. 10. Construction of the “square” [[4, 1,4 ]]q code. Using the 
four physical qubits in the upper and the lower half of the Chimera 
unit cell, we construct two encoded qubits. In (a), the red (thick 
oblique) links represent the penalty couplings (stabilizers), while or¬ 
ange (horizontal), green (vertical) and blue (thin oblique) links form 
the logical Hamiltonian couplings, (b) The section of the encoded 
graph formed by (a). This is also the connectivity induced at level 1, 
where each link represents the collection of level 0 representations in 
(a). 


of the physical qubits comprising each encoded qubit, but we 
use only 2 of these realizations, or alternatively we can use all 
4 but with half the value of h per qubit. Thus 

) = 2he , 45 = 2Jw , (26) 

where we have interpreted the encoding as a level- 1 code. 
This code induces an encoded 2LG graph of level- 1 which 
we can exploit for concatenation. As a second layer of encod¬ 
ing, we use n X n encoded level- 1 qubits to form an encoded 
level-2 qubit on each plane of the 2LG, as depicted in Fig.pT 
Note that with this encoding the 2LG graph structure of level- 
1 is maintained, and that an encoded level-2 qubit comprises 
4n^ physical qubits. 

The level-2 encoding can correct t = 2n^ — 1 errors, i.e., 
has distance A'n? as expected from the code concatenation, and 
will be undecided if there are 2in? errors. To see why this is 
so, consider the level-2 code as a whole, and interpret it as 
a distance d = 4n^ qubit repetition code. A majority vote 
on the level-2 encoded 2LG graph can give the wrong result 
if more than half of the qubits are flipped, i.e., it can recover 
from 2n^ — 1 flips, it is undecided for 2n^ errors, and will give 
the wrong answer for more than 2n^ errors. The level-2 code 
is thus a [[4n^, l,4n^]]o code. 

We note that an additional benefit of the encoding is the 
ability to fine-tune different couplings or local fields in the 
problem Hamiltonian. For example, suppose Jij and hi can 
only take values {— 1,0, !}• How can we generate general 
problems where some couplings have fractional values? To do 
this we use the fact that at the level-2 encoding each coupling 
or local field has 2n different realizations. It is now straight¬ 
forward to rescale all couplings with respect to the largest cou¬ 
pling. E.g., suppose a given coupling J is 9/10 of the largest 



FIG. 11. Level-2 encoding. Green lines represent level-1 couplings 
promoted to penalty couplings (stabilizers), while orange lines rep¬ 
resent level-1 encoded (g) couplings (2 physical couplings). 
Each vertex represents a level-1 encoded qubit (4 physical qubits). 
Red dots represent realizations of level-1 encoded operators. The 
encoding uses n x n level-1 encoded qubits, with the n = 3 case 
depicted here. 


coupling Jmax, then with a sufficiently large n, one can turn 
off 1/10 of the realizations of J and thus achieve the desired 
J/Jmax ratio. The larger n is, the greater the accuracy of this 
ratio. If the quantum annealing device has other possible val¬ 
ues of Jij and hi, then this reduces the necessary value of n 
for a desired ratio. 

It is clear that the square code can be concatenated further. 
The resulting concatenated code can be made to have arbitrar¬ 
ily large distance, i.e., is capable of dealing with increasingly 
more errors, at the cost of more physical qubits per encoded 
qubit. It also exhibits a potentially beneficial growth in the 
encoded local fields and couplings. We describe this in Ap¬ 
pendix]^ 

Decoding of the square code relies on the techniques we 
discussed in Sec. In Appendix we propose an alterna¬ 
tive, recursive decoding scheme. 


VI. EXPERIMENTAL RESULTS 

A. Erustrated-loop Hamiltonians with planted solutions 

We now restrict our attention to a specific set of ran¬ 
dom Ising problem Hamiltonians with “planted” ground states 
lIMl 1821 , which were recently introduced in the context of 
quantum annealing and studied in detail using the DW2 device 
m, albeit without ME and QAC considerations. We briefly 
review how these problem Hamiltonians are constructed; for 
a detailed explanation see Ref. 1341 . 

By “planting” we mean that the planted state is, by con¬ 
struction, a ground state of the problem Hamiltonian. The 
planted state need not be the unique ground state (in fact it is 
typically degenerate), but knowing it means that we also know 
the ground state energy. This advance knowledge circumvents 
the need to resort to computationally expensive exact solvers 
that would otherwise be required to determine such properties 
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of the ground state. Without loss of generality we can always B. Experimental methods 

choose to plant the all-zero state: 


l^-o) = |0i,...,0^) . 


(27) 


We generate random loops (cycles over the Chimera graph) 
with all ferromagnetic couplings, except one randomly cho¬ 
sen anti-ferromagnetic coupling. In this manner each loop is 
a frustrated Ising problem whose ground state is the all-zero 
state of the qubits covered by the loop. A problem instance 
is then built as a sum of such randomly generated frustrated 
loops 0 


h^ = J2hi 


loop 


(28) 


Similarly to the case of constrained satisfiability problems 
(SAT), the hardness of a problem can be tuned by varying 
the number of terms in the sum of Eq. ( [28] ) (i.e., the number 
of clauses in the terminology of SAT). Equivalently, hardness 
depends on a clause density a defined as the ratio between the 
number of loops (clauses) L and the total number of qubits: 


This turns out to generate problems with an easy-hard-easy 
structure as a function of the clause density. We loosely refer 
to the hardness peak as occurring at the critical clause den- 
sity q;c 0 Typically, frustration and hardness peak both at the 
critical clause density; see Appendix[D|for an estimate of frus¬ 
tration. 


Our experiments were performed on the DW2 quantum an¬ 
nealing processor installed at the University of Southern Cal¬ 
ifornia Information Sciences Institute (USC-ISI), which has 
been described in detail before (see, e.g.. Ref. (331). The D- 
Wave devices have been designed in such a way that complete 
logical graphs Kn of n vertices can be minor embedded into 
the Chimera graph representing the physical connectivity [see 
Fig-M |62l|6a. The largest complete graph that can be 
minor embedded into the DW2 hardware graph is K^ 2 - 

Apart from specifying the problem Hamiltonian [Eq. 
the user interface enables us to choose the annealing time 
In all our experiments, the annealing time was set to the small¬ 
est value allowed by the hardware, ta = 20/iS. Eor each 
problem instance we ran 10^ annealing cycles. This was done 
by performing 10 programming cycles of 10^ runs each. Eor 
each programming cycle in which the D-Wave device is pro¬ 
grammed with a given problem Hamiltonian, intrinsic con¬ 
trol errors (ICE) prevent the physical couplings to be set with 
a precision better than ^ 5% of the maximum allowed val¬ 
ues I /imax I , I '/max | • This eiTor iucludes both random (low 
frequency) noise and systematic errors (flux biases). Multi¬ 
ple programming cycles help to average over errors in set¬ 
ting the intended coupling. As an additional precaution de¬ 
signed to average out as much systematic error as possible, 
we performed random gauge transformations between each 
programming cycle. A gauge transformation is an invari¬ 
ance transformation of the Ising Hamiltonian Eq. that is 
achieved by flipping the value of qubit i and by appropriately 
changing the sign of the couplings: 


While the length of each loop is arbitrary, here we only con¬ 
sider loops of length 4 and 6. Note that such short loops were 
not considered in Ref. (^ as they generated instances that 
were too hard (and thus resulted in insufficient solutions to be 
statistically meaningful) with a direct embedding on the full 
(then 503-qubit) Chimera graph. Instead Ref. ED used loops 
of length > 8. Note further that recent work (391 has gener¬ 
ated significantly harder instances for the Chimera graph, and 
we use the current method primarily due to its aforementioned 
convenient feature of advance knowledge of the ground state 
energy. Moreover, while the construction we have described 
obviously generates only a subset of the NP-hard spin glass 
problem instances that can be defined over the 2LG, it fea¬ 
tures the important advantage that the clause density is a tun¬ 
able parameter that allows us to probe the quantum annealer 
on instances of varying and controllable hardness. 


^ In somewhat confusing terminology, this Hamiltonian is “frustration-free” 
in the sense that the ground state of the total Hamiltonian is also the ground 
state of each (frustrated) term in the sum da. 

* We use the term “critical clause density” loosely here, in the sense that 
we do not mean to imply the presence of a phase transition, but merely to 
signal the empirically observed hardness peak. 


o'i —o'i ^ Jij ~Jij hi —hi^ Vj . (30) 


We generated 100 random planted instances per clause den¬ 
sity. Below we plot, for each ct, the mean success proba¬ 
bility and its standard error (one a deviation of the mean) 
for the optimal value of the penalties. Penalties were chosen 
and optimized according to the uniform strategy described in 
Sec. |III A| The penalties were optimized separately for each 
a. Decoding used the strategy of energy minimization over 
the broken qubits (we discuss encoding and decoding strate¬ 
gies in more detail in Sec.jVII Aj). 


We stress that all comparisons between the various strate¬ 
gies (i.e., ME, QAC-ME and D) must properly account for 
the physical resources required by each strategy. Specifically, 
ME requires twice as many physical qubits as D, while the 
level-1 square code (the only case we tested in this work) im¬ 
plementation of QAC-ME requires four times as many. This 
implies that one can run two ME or four D embeddings in par¬ 
allel while consuming the same hardware resources (qubits) as 
needed for a single QAC-ME embedding. To account for this 
parallelism we renormalize and report the experimental D and 
ME ground state (“success”) probabilities as follows: 


Pd 1 - (1 - Pd)^ (31a) 

Pme iH- 1 — (1 — Pme)^ , (31b) 
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(a) Planted (full 2LG graph) (b) Embeddable planted (embeddable 2LG subgraph) 


FIG. 12. Mean success probabilities [after the renormalization of Eq. 0] of instances generated by planting solutions on (a) the square code 
encoded graph, (b) on the embeddable subgraph. Decoding used the strategy of energy minimization. In (a) the comparison shown is between 
the ME of Eig. 1 8(b) I and the QAC-ME of Eig. |8(c)| QAC-ME is seen to improve the ME performance by almost an order of magnitude, with the 
largest improvement obtained at the critical clause density. In (b) the comparison shown is between D of Eig. |9(b)| ME obtained by restricting 
Eig. |8(b)| to the directly embeddable subgraph, and QAC-ME obtained by similarly restricting Eig. |8(c)| It is seen that the performance loss of 
ME with respect to D is completely recovered by QAC-ME. Errors bars denote la confidence intervals here and in all subsequent plots. 


which is the probability of finding the ground state at least 
once after four (D) or two (ME) parallel attempts, using the 
success probability from a single attempt (see Appendixj^for 
a more detailed description of how the quantities above are 
computed from the raw data). 

C. ME vs QAC-ME for the 2LG graph 


compared to D, until a ^ 1.17 is reached. At this value, the 
normalized success probabilities for ME are an order of mag¬ 
nitude smaller than those for the D case. QAC-ME completely 
reverses the performance loss due to ME; the QAC-ME suc¬ 
cess probabilities are comparable to those obtained with the D 
implementation. This again demonstrates that QAC acts as a 
successful error correction strategy. 


Eigure 12(a) shows our experimental results for instances 
with planted solutions as described above and constructed on 
the entire 2LG graph, i.e., the encoded graph of the square 
code shown in Eig. |8(a) Eig. 1 12(a)] compares the ME [cor¬ 
responding to Eig. 8(b)| and QAC-ME [corresponding to 


Eig. 8(c)| implementations on the same set of instances. It 
shows that QAC-ME is effective in boosting the performance 
of the quantum annealer on the entire range of clause densi¬ 
ties. Interestingly, the harder the problems, the larger is the 
advantage of QAC-ME over ME. We empirically identify the 
critical clause density as acrit ~ 0.94, where the success 
probability has a minimum. At this value QAC-ME yields 
success probabilities that are about one order of magnitude 
higher than ME’s. Thus error correction becomes more effec¬ 
tive for harder problems, as was also observed in the context 
of random Ising problems in Ref. IbOl . 


D. D vs ME vs QAC-ME 

We next compare between D and ME by constructing in¬ 
stances with planted solutions on the 2LG subgraph shown in 
Eig. |9(a)[ Experimental results are shown in Eig. 1 12(b)] As ex¬ 
pected, the ME strategy suffers a steady decay in performance 


We note that there is no clear evidence of a critical clause 
density (or hardness peak) in this case. The ME strategy ex¬ 
hibits a plateau after a ^ 1.17, while the D and QAC-ME 
success probabilities continue to decrease, with QAC-ME de¬ 
creasing more slowly. Up to this value, the increase in hard¬ 
ness is intrinsic and can be related to the increase in frus¬ 
tration. It may be possible to associate the further increase 
in hardness for larger clause densities to a precision effect, 
e.g., an increase in the range (the ratio between the largest 
and smallest couplings in the problem Hamiltonians; see Ap¬ 
pendix]^ for a quantification of this effect). This adversely 
affects the performance of the DW2 because a large coupling 
range means that some couplings are set to smaller values 
than can be physically implemented (the DW2 has a range 
of about 4 bits), due to the well-studied effect of ICE that 
leads to the specification of the wrong problem Hamiltonian 
O [331 ESI H [73. The saturation of the ME success prob¬ 
abilities can be understood as being due to having arrived at 
this precision limit before D and QAC-ME, i.e., for a > 1.17 
the ME results are dominated by ICE, and the improvement 
due to QAC is at least in part due to countering ICE. The D 
implementation is, of course, also susceptible to ICE, but is 
impacted less than ME due to the absence of broken qubits in 
D. 
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FIG. 13. Success probabilities of ME and QAC-ME for uniform and nonuniform penalties for the four types of problem instances, normalized 
to the success probability of ME with uniform penalties. Decoding used the strategy of energy minimization. In all cases ME with uniform 
penalties performs worst, across the entire range of clause densities. ME with nonuniform penalties improves the ME performance somewhat, 
but always less than QAC-ME with uniform penalties. QAC-ME with nonuniform penalties does worse than QAC-ME with uniform penalties 
when the instances have uniform couplings [(a),(c)], but almost always does better than QAC-ME with uniform penalties when the instances 
have non-uniform couplings, as in (b) and (d). In the deformed embeddable case (d) QAC-ME with nonuniform couplings does significantly 
better even than the direct embedding. 


E. Optimized uniform vs nonuniform penalties 

In Sec.|I^we discussed optimizing the energy penalties and 
proposed a new nonuniform strategy. In this section we exper¬ 
imentally test the effectiveness of this strategy and compare 
it to the previously proposed uniform strategy. Specifically, 
we compare the performance of the square code when penal¬ 
ties are chosen according to either Eq. (Tf) (uniform case) or 
Eq.® (nonuniform case), where the overall penalty energy 
scale 7 is a parameter that must be optimized independently. 

Intuitively, a nonuniform choice of penalties should per¬ 
form best when the logical couplings are distributed nonuni- 
formly on the logical graph. The problem instances we have 
considered so far were generated by adding frustrated loops 
that are uniformly distributed over the logical graph, resulting 


in similarly uniformly distributed couplings. We thus consider 
two additional sets of nonuniform instances. 


1. Weighted planted frustrated-loop instances 

The first set of instances is a simple modification of the con¬ 
struction using frustration-free Hamiltonians described above. 
The only difference compared to Eq. ( [^ is that here we con¬ 
sider a weighted sum of the same set of frustrated loops as 
used to generate Eig. with the weights varying nonuni- 
formly on the logical graph: 

L 

= ( 32 ) 
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(a) Planted, uniform 


(b) Planted, nonuniform 




(c) Weighted planted, uniform 


(d) Weighted planted, nonuniform 


FIG. 14. Comparison between different decoding strategies for the ME and QAC-ME implementations, for planted [(a) and (b)] and weighted 
planted [(c) and (d)] instances. Panel (a) and (c) [(b) and (d)] is for the case of uniform [nonuniform] penalties. The success probabilities plotted 
are normalized to the ME case with random decoding (Pmect)- Panels demonstrate the improvement of the decoded results when going from 
CT to MV and EM. In the ME case, the EM strategy resoundingly beats CT. In the QAC-ME case the ordering is EM>MV-EM>MV-CT>CT, 
with the advantage of the best strategy (EM) being more significant in the nonuniform case. 


We choose the Wi as follows. We define yi^Zi} to be the 
“center” of loop i on the 2LG (defined as the integer part of 
the mean of the coordinates of the spins comprising the loop), 
and set Wi = Xi. This increases the typical strength of the 
couplings by a factor of 8 between the two opposite sides of 
the 8x8x2 2LG logical graph embedded into the DW2. This 
introduces a spatial non-uniformity and also reduces the mag¬ 
nitude of the smallest couplings of this set of instances (see 
Appendix]^, thus making the instances significantly harder, 
as we shall see below. Note that the same planted solution 
as before (|^o)) applies, independently of the choice of the 
weights. 


2. Deformed embeddable instances 


The second set are nonuniform instances we generate by 
“deformation” of the set of instances with planted solutions 
defined on the embeddable subgraph [Fig. 9(a)|, as follows. 
For each instance we randomly pick a subset / of encoded 
qubits (vertices of G) and apply a linear transformation (shift¬ 
ing and compression) to the encoded couplings such that: 


Ji,j yi E ^ I (33a) 

J-j ^ (2sign(J-j) + J-j)/3 Vi, J ^ 7 . (33b) 


This partitions the encoded qubits into three sets 7, J, 
where the set 7 (the one picked) is connected only through 
coupling that have been transformed according to Eq. ( |33a| ) 
(small couplings), K is connected through couplings that have 
been transformed according to Eq. ( |33b| ) (large couplings), 
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and J is connected through couplings that have been trans¬ 
formed according to both Eqs. ( |33a| ) and ( |33b| ). 

The average of the encoded couplings is ^ 0.5. Using 
this value in Eqs. ( |33a| ) and ( |33b| ) yields average coupling 
strengths of 1 /6,1/2 and 5/6 for the sets 7, J, K respectively. 
We randomly pick 35 qubits out of the available 120 for the 
I set, which yields three similarly sized sets. Under this de¬ 
formation, the originally planted configuration is no longer 
necessarily the ground state. In this case, we use the exact 
solver that uses a bucket tree elimination strategy (BTE) and 
is included in the D-Wave user libraries to determine the exact 
ground state energy. This exact solver is specifically designed 
to solve small. Chimera-structured problems. 

Eigure[^ shows the experimental results for ME and QAC- 
ME, with both uniform and nonuniform penalties. The panels 
show the four sets of instances: planted and weighted planted 
(both defined on the full 2LG), embeddable planted and de¬ 
formed embeddable (defined on the embeddable subgraph). 
In all panels of Eig. we see that QAC-ME outperforms 
ME. The choice of nonuniform penalties improves the suc¬ 
cess probabilities of both the ME and QAC-ME in the case 
of nonuniform instances [panels (b) and (d)]. Moreover, we 
see that the nonuniform choice of penalties improves the ME 
implementation in the uniform case too, while it does not im¬ 
prove the QAC-ME scheme, suggesting that QAC-ME with 
a uniform penalty is already close to optimal in the case of 
uniform instances. (The dependence of the optimal penalty 
strength on clause density is shown in Appendix]^) 

As a final remark of this section, note the large boost of 
nearly two orders in magnitude in the success probabilities 
that QAC-ME achieves over ME in the weighted case shown 
in Eig. |13(b)[ In this case QAC-ME is particularly success¬ 
ful. A similar advantage of QAC-ME over ME is observed 
in Eig. 1 13(d)] where the advantage of QAC-ME with nonuni¬ 
form penalties is the most pronounced, easily outperforming 
the direct embedding as well [in contrast to the tie in the em¬ 
beddable planted case, seen in Eig. 12(b)| . 


VII. IMPACT OF THE DECODING STRATEGY 


We now consider the impact of various decoding strate¬ 


gies defined in Sec. IV on performance as measured in terms 
of success probabilities. As already discussed in Sec. |IVB 


because both the ME and QAC-ME implementations of the 
square code involve an even number of physical qubits (two 
and four respectively), the possibility of ties means that bro¬ 
ken qubits can be undecodable by majority vote. In the ME 
case, in particular, every broken qubit is indeed a tie. In the 
QAC-ME case there are two possibilities: “partially” broken 
qubits (decodable through majority vote) and ties. 

We thus consider the following decoding strategies: 


1. ME: 


(a) Decode ties randomly via coin tossing (CT); 

(b) Decode ties via energy minimization (EM). 

2. QAC-ME: 


(a) Decode all broken qubits via CT; 

(b) Decode all broken qubits via energy minimization 
(EM); 

(c) Decode partially broken qubits via MV and ties 
via CT (MV-CT); 

(d) Decode partially broken qubits via MV and ties 
via EM (MV-EM); 

(e) Decode partially broken qubits via MV and ties 
via EM. Randomly select a number of unbroken 
qubits equal to the number of partially broken 
qubits and decode via EM (MV-EM(R)). 


A. Performance of the various decoding strategies 


While we do not expect random decoding (the CT strat¬ 
egy) to work well, it sets a useful baseline against which we 
can compare our other decoding strategies. The comparison 
between all strategies in the case of instances defined on the 
full 2LG graph is shown in Pigs. and Pigure[^is for 
the planted and weighted planted instances and reports results 
for both the ME and QAC-ME cases. Pigure focuses on 
the MV-EM(R) decoding strategy (QAC-ME only). Comple¬ 
mentary results for the embeddable planted and the deformed 
embeddable cases are shown in Appendix [G] 


As expected, CT is the least effective decoding strategy, and 
is substantially outperformed by MV. The optimal choice is 
always the EM strategy (recall Sec. |IV C| ), which yields im¬ 
provements over both MV and MV-EM. The benefit of EM 
is larger in the case of nonuniform penalties and weighted in¬ 
stances. In Pig. 14(d)| we see that EM decoding improves the 
success probabilities by a factor of about 3 as compared to the 
second best decoding strategy, MV-EM. The advantage of EM 
over MV and CT is directly related to the number of broken 
qubits at the optimal penalty value. The larger this number, 
the more bit flip errors we expect to correct during the decod¬ 
ing process by doing energy minimization. As we shall see in 
the next section, the typical number of broken qubits is usually 
larger for the nonuniform penalties and for the nonuniform in¬ 
stances we have constructed. 


It is interesting to check whether unbroken qubits (where 
all physical qubits agree), which cannot be decoded via MV, 
can be further improved by EM. The MV-EM(R) decoding 
strategy is designed to probe this, by singling out a fraction 
of the unbroken qubits and applying the EM strategy to them. 
We thus compare MV-EM(R) to MV-EM and EM decoding 
in Pig. 15 again for the planted and weighted planted in¬ 
stances defined on the full 2LG graph, with nonuniform penal¬ 
ties. We find that the MV-EM(R) strategy offers no advantage 
over MV-EM, leaving EM as the superior strategy. Since MV- 
EM(R) is essentially indistinguishable from MV-EM, we may 
conclude that the unbroken qubits make a negligible contribu¬ 
tion to the error rate. 
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(a) Planted, nonuniform (b) Weighted planted, nonuniform 


FIG. 15. Comparison between different decoding strategies for the QAC-ME implementation with nonuniform penalties, for (a) planted and 
(b) weighted planted instances. The MV-EM(R) strategy, where a random set of unbroken qubits is subjected to energy minimization, is 
shown along with MV-EM and EM. It can be seen that MV-EM(R) does no better than MV-EM. This implies that logical errors are mostly 
concentrated on broken qubits. The EM data shown here is the same as for the QAC-ME case in Eigs. |14(b')] and |14(d^ 
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EIG. 16. Broken qubit median probabilities for the “planted” in¬ 
stance set at the critical clause density a = 0.94. Results are shown 
for the full set of 100 instances implemented with the QAC-ME 
scheme, for the optimal value of the (uniform) penalty 7 = 0.2. The 
abscissa is the index of the encoded qubits on the 2LG graph shown 
in Eig.| 8 (a)| 


B. Broken qubit probability and the percolation threshold for 
efficient decoding 

Following the discussion in Sec. |IVD] efficient optimal de¬ 
coding via the EM strategy is possible if the probability pbq 
to find a broken encoded qubit is smaller than the percolation 
threshold P 2 LG of the 2LG. This ensures that no large clusters 
of broken qubits are formed that would make energy mini¬ 
mization too computationally costly. While we do not know 
the value of the percolation threshold of the (infinite) 2LG, it 
must certainly be bounded between the thresholds Pcubic and 


Psquare fof cubic and Square lattices, respectively |[84ll85l : 

0.3116 = Pcubic PlLG ^ Psquare — 0.5927 . (34) 

We may further expect P 2 lg to be closer to Psquare than to Pcubic 
since the 2LG is geometrically more similar to a square lattice. 
Recall that the 2LG is the encoded graph of the square code, 
and hence we identify the latter’s percolation threshold with 

P2LG- 

We first check that the experimental probability pbq to find 
broken qubits in the square code depends only weakly on the 
position of the encoded qubits. This is confirmed (for the 
QAC-ME case) in Eig. at the optimal penalty value for a 
fixed clause density, and justifies the assumption underlying 
our percolation threshold argument about the efficiency of the 
EM decoding strategy in Sec. |IV D| 

Eigure[^ shows the probability of a broken qubit as a func¬ 
tion of the clause density, for both ME and QAC-ME, in the 
case of instances with planted and weighted planted solutions 
on the full 2LG, for both uniform and nonuniform penalties. 
The most important finding is that in all cases, the experi¬ 
mental probability of a broken logical or encoded qubit is 
well below Pcubic. and hence below the percolation threshold 
of the square code. This suggests that scalable and efficient 
decoding of the square code through energy minimization is 
possible, at least for the (frustrated) class of instances con¬ 
sidered here. Additional supporting arguments are given in 
Appendix [H| using a worst-case analysis. 

There are a few other interesting observations we can make 
about Eig.[^ Eirst, the very fact that the observed pbq is non¬ 
vanishing is the reason that a decoding of ME and QAC-ME 
is necessary. Second, note that in both the ME and QAC-ME 
cases, the optimal penalty values tend to keep the probability 
Pbq of broken qubits fairly constant over the range of clause 
densities a we have studied. The optimal penalty values can 
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FIG. 17. Probability pbq of broken qubits in ME and QAC-ME, as a function of the clause density, for planted (top) and weighted planted 
(bottom) problem instances defined on the full 2LG graph. Left: uniform penalties. Right: nonuniform penalties. ME almost always has a 
lower pbq than QAC-ME. Also shown is the probability of tied qubits in the QAC-ME case (which is already included in the broken qubits 
probability). 


thus be converted into an optimal broken qubits probability 
that depends on the set of instances chosen and on the im¬ 
plementation (ME or QAC-ME) considered. Third, note that 
Pbq is typically smaller for ME than for QAC-ME, which 
is unsurprising given that the QAC step in QAC-ME adds a 
second way for qubits to break, in addition to ME. Eourth, 
note that pbq is typically larger with nonuniform penalties 
or weighted instances. The intuitive reason for this is that 
nonuniform penalties are by design weaker for qubits that are 
less strongly coupled. This facilitates the appearance of more 
broken qubits. This is consistent with the experimental find¬ 
ing of the previous section that EM decoding results in a larger 
advantage in these cases, simply because there are more bro¬ 
ken qubits to be optimized. Einally, note that the differences 
in performance of the various decoding strategies considered 
in Sec. VIIA can be related to pbq- 


VIII. TESTING TEMPERATURE AND NOISE 
DEPENDENCE VIA SIMULATED QUANTUM ANNEALING 


To gain a better understanding of the mechanisms that are 
most responsible for the experimental success of QAC, we 
conducted numerical simulations that illuminate the role of 
temperature and ICE. A master equation simulation treating 
the device as an open quantum system 1^ is unfeasible due 
to the large number of qubits involved in our experiments. 
Hence we resort to a semi-dynamical approach based on quan¬ 
tum Monte Carlo. The approach, called simulated quantum 
annealing (SQA) (Tt] [86l [87l, attempts to mimic the evolu¬ 
tion of an open quantum system where stochastic processes 
(typical of a noisy regime) dominate the dynamics, and es¬ 
sentially provides the instantaneous Gibbs state along the an¬ 
nealing evolution (see Appendix |I] for details). This approach 
has successfully reproduced the success probability distribu¬ 
tion of the D-Wave devices for random Ising instances (TtI, 
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(a) T = 2, X = 5% (b) T = 2, X = 10% 



(c) T = 4, X = 5% 


FIG. 18. SQA simulation results for the “embeddable planted” instances class, (a) Relatively low temperature T and low noise strength x^ 
QAC-ME does not provide an advantage over the uncorrected ME approach. Note how success probabilities are typically larger for SQA in 
this case with respect to the experimental results [Eig. | 12(b)) . (b): Relatively low temperature and high noise strength: QAC-ME gains a clear 
advantage over ME but is still inferior to D. (c) Relatively high temperature and low noise strength: QAC-ME gives a consistent advantage 
over ME, and it is competitive with D. This situation is the most similar to the experimental results. We used 20,000 sweeps in all cases shown. 
Temperature is in units of B{tf)/kB- 


although it does not provide a complete description when one 
also accounts for ground state degeneracy and excited states 
ED, or open system dynamics . 

Our strategy is not to use SQA as a faithful model for the 
DW2 device. We thus do not attempt to tune the SQA param¬ 
eters to find a best fit for the experimental results. Rather, we 
use SQA as a numerical tool to explore the response of our 
QAC-ME scheme to changes of various important quantities 
that cannot be controlled experimentally, specifically temper¬ 
ature and noise. The noise level x is the standard deviation of 
statistically independent Gaussian noise added to all physical 
couplings Jij and local fields hi, and is introduced to repro¬ 
duce ICE. We performed simulations for different values of 
the number of sweeps (a sweep is a single update of all spins), 
but found that the results have a very weak dependence on this 
number, so it is not a crucial parameter in our analysis. On 


the other hand, the temperature and the noise level heavily af¬ 
fect our numerical results. Eiguresummarizes our numer¬ 
ical results for the “embeddable planted” instances [compare 
with the experimental results shown in Eig. 12(b)| . It shows 
that QAC-ME becomes competitive with respect to ME and 
D only for sufficiently large temperature and/or noise. We re¬ 
gard this result as strong evidence that the QAC schemes are 
particularly efficient in reducing the amount of errors due to 
thermal noise and imprecisions in parameter settings. 


IX. CONCLUSIONS AND DISCUSSION 


In order to exploit the potential advantages of quantum an¬ 
nealing over classical solvers in solving “non-native” opti¬ 
mization problems they must accommodate a minor embed- 
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ding (ME) step that maps the given problem to the given 
device hardware graph. In this work we introduced a new 
quantum annealing correction (QAC) scheme that is compati¬ 
ble with the minor embedding step—the “square code”—and 
studied its efficacy experimentally using the DW2 device. We 
demonstrated that this hybrid QAC-ME scheme significantly 
boosts performance by applying it to frustrated Ising problems 
with known (“planted”) ground states. In the construction 
of such instances, a controllable clause density parameter a 
tunes both hardness and frustration. The largest performance 
boost (measured in terms of the success probability of finding 
the ground state) via QAC-ME occurs at the critical clause 
density where hardness peaks. 

We have also considered problem instances that are not 
only embeddable into the square code but also directly into 
the physical “Chimera” hardware graph of the DW2 device. 
This allows for a comparison of the direct, ME and QAC-ME 
embeddings for the same problem instance. As expected ME 
behaves poorly with respect to the direct embedding. The per¬ 
formance improvement achieved with QAC-ME is particular 
significant at o^c and suffices to match that of the direct em¬ 
bedding. This is a proof-of-concept that QAC-ME can (and 
should) be used to overcome the limitations of minor embed¬ 
ding on quantum annealers. QAC-ME thus has the potential 
to extend the competitiveness of these devices beyond natively 
embeddable optimization problems. 

Several questions remain to be addressed in future work, 
requiring the use of next-generations quantum annealers. The 
most important of these is arguably the scaling of error- 
corrected quantum annealing. To properly address this ques¬ 
tion we need to be able to embed larger logical problems (the 
square code considered in this work implements up to 120 
encoded qubits) on devices with a larger number of physical 
qubits. An important aspect for any scaling study is the deter¬ 
mination of the “optimal” implementation. In the case of a di¬ 
rect embedding, this boils down to the question of determining 
the optimal annealing time (see, e.g.. Ref. |[34l ). Optimality 
in minor embedded quantum annealing is a broader question. 
This is due to the freedom in choosing the minor embedding 
and, once this is given, the exact map between the logical and 
the physical couplings. In particular, the performance of mi¬ 
nor embedded quantum annealing is sensitive to the choice of 
energy penalties, and determining the optimal choice of penal¬ 
ties is important both for determining the scaling of physi¬ 
cal quantum annealing and for optimizing the performance on 
practical applications. We have shown that choosing nonuni¬ 
form penalties that depend on the local strength of the logical 
couplings can give a significant advantage with respect to a 
uniform choice. 

The optimal choice of the decoding procedure is another 
important question. We found global energy minimization 
over broken qubits to be the optimal choice for our problem 
instances. We showed that in order for this decoding tech¬ 
nique to be performed at a negligible computational cost, the 


typical probability to have broken logical or encoded qubits 
should be lower than the per-site percolation threshold of the 
logical or encoded graph. It turns out that the square code, for 
the frustrated problems we have considered, is efficiently and 
optimally decodable through energy minimization even when 
its size is scaled. Our results demonstrate that the square code 
on the tested DW2 lies in a “decodable” regime, with a suf¬ 
ficiently small number of broken qubits when the penalties 
are optimized. It is important to study whether optimal de¬ 
coding can be performed efficiently in other ME and QAC- 
ME schemes. Whenever this is not possible, efficient, albeit 
suboptimal, strategies should be developed. To this end, it 
is likely that techniques developed in the decoding of Low 
Density Parity Check (LDPC) codes |[89l will find important 
applications in quantum annealing too. In particular, subop¬ 
timal decoding through belief propagation 1901 is a promis¬ 
ing approach to successfully decode ME quantum annealing 
when the broken qubit probability is larger than the decoding 
threshold. 

More generally, it is important to understand how the per¬ 
formance of minor embedded quantum annealing degrades 
with the increase of the size of the logical or encoded qubits. It 
is likely that the larger the logical qubits required for a given 
minor embedding are (this will be the case for increasingly 
larger logical problems) the more important QAC-ME will 
be as a tool to prevent the degradation of the performance of 
quantum annealers on minor embedded implementations. 

While our experimental conclusions are naturally limited 
to the particular set of problem instances we studied, we may 
extrapolate that our results more generally support the effec¬ 
tiveness of the QAC-ME strategy. Confirming this can have 
important implications in scaling quantum annealing to solve 
larger and harder optimization problems of practical interest, 
and we expect QAC-ME techniques to play a central role in 
achieving this goal. 
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Appendix A: Energy Minimization Decoding 


level-r fields and couplings are 


We performed decoding through energy minimization by 
implementing simulated annealing on the “decoding” Hamil¬ 
tonian defined in Eq. ( |2Q| ). We used a linear schedule for the 
temperature, with the initial and final temperatures equal to 
Tin = 4max^j(J^j) and T^n = 0.1 min^j( J^j) respectively. 
We chose the computational effort in the decoding step to be 
constant for all the decoding attempts. For each readout after 
an annealing cycle, decoding consisted of repeating 10 simu¬ 
lated annealing runs with 10 sweeps. For each decoding at¬ 
tempt, we defined a decoding success probability Pdec as the 
fraction of times the lowest decoded configuration was found. 
We consider this quantity as a measure of our decoding qual¬ 
ity, with a Pdec close to 1 meaning that the best decoded con¬ 
figuration is likely to be the optimally decoded state. A small 
value of Pdec 5 on the other hand, means that a better decoded 
configuration could likely be found with increased computa¬ 
tional effort (a larger number of sweeps and/or repetitions in 
the simulated annealing procedure). 

Pdec is shown in Fig.[^for a few representative examples. 
Recall that the number of broken qubits is typically larger for 
QAC-ME (Fig. [^, which therefore requires a larger decod¬ 
ing effort in order to obtain a similar decoding quality. This 
explains why we observe lower decoding success probabilities 
for the QAC-ME case with energy minimization. Nonuniform 
penalties are also responsible for a larger number of broken 
qubits. Improving Pdec for QAC-ME by expending more ef¬ 
fort on the simulated annealing step would thus further in¬ 
crease the advantage of QAC-ME with respect to ME, as well 
as the advantages of nonuniform penalties with respect to the 
uniform case. 


Appendix B: Concatenating the square code 


In Sec. V C of the main text we mentioned that the square 
code can be concatenated, and described the first and sec¬ 
ond encoding levels. To see how this can be continued, con¬ 
sider, e.g., combining n x n encoded level-2 qubits to form 
an encoded level-3 qubit on each plane of the 2LG, to ob¬ 
tain a [[4n^, 1, 4n^]]o code. In general we have [n,r] = 
1,for a level-r code with integer r > 
1. The realizable stabilizer terms included in the penalty 
Hamiltonian suffice to penalize every possible bit-flip error, 
at every level of concatenation. 

It is interesting to consider the growth in the encoded fields 
and couplings. As illustrated in Fig. [m for levels r > 2 
the encoded couplings grow in proportion to the side length 
of the square of encoded qubits, and the encoded fields must 
grow in the same proportion in order to maintain their rela¬ 
tive strengths. E.g., the level-2 local fields correspond to n 
level-1 realizations, each with two level-0 realizations. Anal¬ 
ogously, every encoded level-2 coupling comprises n level-1 
couplings, which are in turn equivalent to 2 level-0 couplings. 
Thus = 2nh£ and = nJ^l) = 2nJi£>. 

At level r each square has side length so the encoded 


= 2n^-^ht , J'fj = 2vJ-^Ju' , (Bl) 

i.e., the energy scale of the level-r final Hamiltonian has been 
boosted by a factor of 2rf~^. Note that the encoded penalty 
strength grows in proportion to the area of each square if all 
the penalties are utilized. Since this would result in overly 
rigid qubits a simple remedy is to use only a linear number of 
the total available penalties. 

Since the number of physical qubits required grows 
quadratically faster than the strength of the encoded local 
fields and couplings, it is not a priori clear that a higher con¬ 
catenation level r is beneficial. However, in a Markovian 
model where the probability of a thermal excitation at inverse 
temperature p is proportional to the Gibbs factor (e.g., 
Ref. 1^ ). with A the gap from the ground state, it is possible 
to realize a gain by increasing r, provided the gap grows along 
with the energy scale of the level-r final Hamiltonian. Deter¬ 
mining the scaling of the gap with r is a non-trivial problem 
since the energy scale of the initial transverse field Hamilto¬ 
nian does not also grow under the present encoding, as re¬ 
marked earlier. Our proposal for concatenating the square 
code for QAC can be viewed as a physically realizable (and 
hence severely compromised) version of Bacon’s proposal for 
quantum concatenated-code Hamiltonians ED. 


Appendix C: Recursive decoding of the square code 

In this section we propose an alternative (untested) recur¬ 
sive decoding scheme of the square code, that is similar to the 
energy minimization strategy. 

Suppose we already decoded all the untied encoded qubits 
with only the tied encoded qubits remaining to be decoded. 
Consider a problem Hamiltonian j)eG JijCrfcTj, 

where the sum is over all connected pairs of logical qubits i 
and j in the encoded graph G (the graph obtained after, say, 
the first level of encoding using the square code), and for no- 
tational simplicity here we use undecorated Pauli operators 
to denote the encoded operators. At the end of a (quantum) 
annealing run some encoded qubits will be tied. The untied 
qubits now have fixed values {s^}, so they act as local fields 
on the tied qubits, weighted by the couplings. We can rewrite 
the problem Hamiltonian after this run as 

= T + E > (Cl) 

j'eties 

where the sum runs only over the tied qubits, which occupy 
the vertices of a subgraph G' of G, and where the local fields 
are 


hj' — ^ ^ Jij'Si , (C2) 

i 

with the sum including all the untied qubits coupled to the 
given tied qubit j'. 
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The renormalized Hamiltonian represents a new Ising 
problem. To find the values of the tied qubits we can run 
on the (quantum) annealer, i.e., implement the EM 
strategy by bootstrapping the problem. If after such a run 
there are still tied qubits we can repeat the same process, 
defining a new Hamiltonian + 

with hj» = etc. 

Note that successive runs of the quantum annealer need not 
necessarily resolve the tied qubits, but since G" is by con¬ 
struction a subgraph of G', this process will converge. Any 
ties that remain after the process has converged can be de¬ 
coded using the previously described EM strategy. Account¬ 
ing for wall-clock times, it may be preferable to start running 
this sequence of recursive optimizations on a classical solver 
once the subgraph size has become small enough. 


Appendix D: Broken Qubit Probability, Range and Frustration 


In this section we present the complete data set in Eigs.[20| 
(planted), (embeddable planted), (weighted planted), 
and [^(deformed embeddable) for the broken qubit (and ties) 
probabilities for the instances studied for all values of the 
penalties, for both the ME and QAC-ME cases. Note how, as 
expected, the broken qubit probability (pbq) strongly depends 
on the strength of the energy penalties. One more feature that 
is evident from the same figures is the different behavior of 
Pbq in the uniform and nonuniform cases, at a given value 
of the penalty strength 7. In the uniform case, pbq is mostly 
a monotonically decreasing function of the clause density a, 
while in the nonuniform case there is a clear peak that is lo¬ 
cated close to the critical clause density. 

There is a simple and intuitive explanation of this behavior 
in terms of range and frustration, that is consistent with the 
general discussion of Sec. IIP Let us first give a quantitative 
definition of both frustration and range. We define frustration 
as the percentage of Jij couplings that are not satisfied on 
the planted ground state solution. The inverse range is the 
smallest (non-zero) absolute value of the Jij couplings of a 
given instance. Each time we implement a given instance, 
we always set the value of the largest coupling equal to 1. 
This exploits the largest energy scale achievable with the DW2 
processor. The range thus gives us an estimate of the required 
precision needed to correctly implement the given instance. 

Eigure shows the computed values of the required (in¬ 
verse) range and the frustration of the various classes of in¬ 
stances. We can see that the qualitative behavior of pbq in the 
uniform case follows that of the inverse range. This is related 
to the fact that a smaller inverse range implies the presence 
of a larger number of weak couplings that are dominated by 
the uniform and relatively larger energy penalty, making the 
appearance of a broken qubit less likely. On the other hand, 
the value of the energy penalties in the nonuniform case is ad¬ 
justed to the typical value of the logical couplings. Typically, 
a smaller inverse range results in weaker energy penalties, bal¬ 
ancing the range effect. Thus, in the nonuniform case pbq is 
more positively correlated with frustration. This is consistent 
with the general expectation that pbq is large for highly frus¬ 


trated instances. 


Appendix E: Data Collection and Analysis 


As already mentioned in Section [ylB] the presence of con¬ 
trol errors (ICE) due to low frequency noise and systematic bi¬ 
ases prevents the physical couplings to be set with a precision 
better than ^ 5%. To average out the effects of such errors, we 
run 10 programming cycles for all the encoding schemes (D, 
ME, QAC-ME) and penalty values for every instance consid¬ 
ered. Each programming cycle implements a gauge transfor¬ 
mation and consists of 10^ annealing runs. Eor each instance, 
encoding strategy and penalty value, (altogether labeled here 
by the index i) we extract an “intrinsic” success probability P* 
as the overall probability to find the ground state over the 10 
programming cycles. In other words, this is the mean of the 
success probabilities extracted from each programming cycle, 
and it serves as our best guess for the intrinsic (without ICE) 
success probability of each encoding of a given instance 






(El) 


We then use these gauge-averaged probabilities to compute 
the renormalized success probabilities for the D and and ME 
cases: 


^ 1 - (1 - , (E2a) 

^ 1 - (1 - , (E2b) 

Eqac-me Pqac-me ■ (E2c) 

As explained in the main text, this renormalization is neces¬ 
sary to correctly take into account the amount of hardware re¬ 
sources that each implementation requires. We then use these 
renormalized probabilities to compute 5000 means, obtained 
from 5000 bootstrapped samples of 100 instances (picked 
with replacement among the 100 total number of instances) 
for each clause density a and penalty value 7. The quantities 
plotted in Eig.[^(as well as in all other figures plotting ex¬ 
perimental data) are the mean and the standard error of such 
5000 samples of bootstrapped means. 

The same approach was followed to compute means and 
error bars of other similar plots. See, e.g., the plot showing 
broken qubits probabilities (Eig.p^. 


Appendix F: Energy Penalties 

We show explicitly in Eig.j^the value of the optimal 7 for 
all the cases considered in the main text. As can be seen, this 
value is fairly constant as a function of the clause density a. 
There is a slight tendency for the optimal 7 to decrease with 
a in the case of uniform penalties. This is most pronounced 
in the ME implementation of the deformed embeddable in¬ 
stances. When choosing nonuniform penalties, however, the 
optimal 7 again becomes almost constant. This shows that 
the nonuniform choice of penalties has the potential of being 
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closer to optimality than the uniform case with a single choice 
of the optimal penalty strength. Since optimization of the en¬ 
ergy penalties can be time consuming, this property is impor¬ 
tant for any practical encoding scheme. Ideally, we would like 
an encoding scheme (for the choice of penalties) that is close 
to optimality and does not depend on any quantity that has to 
be optimized on an instance-by-instance case. 


Appendix G: Decoding Strategies 


We show in Figs. 26 and 27 the comparison between var¬ 
ious decoding strategies for the embeddable planted and the 
deformed embeddable cases. This complements the results 
shown in Figs. and [T^ of the main text, which showed the 
results for the planted and weighted planted instances. Results 
for the embeddable sets are very similar to the sets defined on 
the full 2LG. We notice here too that the use of the EM strat¬ 
egy is crucial for the QAC-ME encoding to give comparable 
(for the planted embeddable set) or better (for the deformed 
embeddable set) results than the direct implementation (D). 


Appendix H: Worst-case Scenario Decodability of the Square 
Code 

We can set an upper limit (worst-case scenario) for the 
probability to obtain ties and broken qubits in general by as¬ 
suming that all four physical qubits comprising an encoded 
qubit are uncorrelated. This corresponds to assuming that the 
output of the annealer corresponds to a thermal distribution 
with infinite temperature. Within a single encoded qubit there 
are 6 physical qubit states that correspond to a tie (out of a 
total of 16 states), so we obtain = 6/16 = 0.375 for the 
worst case probability to obtain a tie. The worst case probabil¬ 
ity to obtain a broken qubit is = 14/16 = 0.875. While 
p^^ is definitely above the percolation threshold of the square 
code, p^l is almost certainly below it since the device out¬ 
put is better approximated by a finite temperature thermal dis¬ 
tribution (as opposed to an infinite temperature distribution). 
Moreover, the non-vanishing penalty term and (if present) lon¬ 
gitudinal field disfavor broken qubits and ties. This strongly 
suggests that in the more realistic scenario the number of bro¬ 
ken qubits is well below the worst-case value. Fig. shows 
the mean size of the connected domains over 10^ randomly 
generated configurations for increasingly large N x N 2LG 
graphs. The fact that this mean size saturates instead of grow¬ 
ing with N is evidence that p^l < P 2 LG- This means that the 
square code is efficiently decodable (at least using using the 
MV-EM strategy) even in the worst-case scenario. 


Appendix I: Simulated Quantum Annealing 


We reported SQA results in Sec. VIII of the main text. Here 
we briefiy review this technique. SQA is a quantum Monte 
Carlo based algorithm whereby Monte Carlo dynamics are 
used to sample from the instantaneous Gibbs state associated 
with the Hamiltonian in Eq. Q- The state at the end of the 
quantum Monte Carlo simulation of the quantum Hamilto¬ 
nian H(t) is introduced as the initial state of the Monte Carlo 
simulation with Hamiltonian H{t At). This proceeds un¬ 
til H(tf) is reached. Originally proposed as an optimization 
algorithm |[86l [87l, it has gained some traction as a classi¬ 
cally efficient model for finite temperature quantum annealers 
ifTTl [2n [33l [34l. SQA does not capture the unitary dynam¬ 
ics of the quantum system, but it is hoped that the sampling 
of the instantaneous Gibbs state captures thermal processes in 
the quantum annealer, which may be the dominant dynamics 
if the evolution is sufficiently slow. Although there is strong 
evidence that SQA does not completely capture the final-time 
statistics of the D-Wave processors |[2T]|28l, at present it is the 
only viable means to simulate large open quantum annealing 
systems. 


We used discrete-time quantum Monte Carlo in our simula¬ 
tions, which we briefiy review. For Hamiltonians of the form 
of the transverse Ising model, such as in Eqs. Q and ([^, the 
sampling from the instantaneous Gibbs state is done by sam¬ 
pling from the dual classical system with Hamiltonian He and 
inverse temperature P : 

E (II) 

i,r 

i i<j 

where Nr is the number of Trotter slices used along the Trotter 
direction, Si^r denotes the ith classical spin on the rth Trotter 
slice, and 




J±{t) = - In [tanh(/3A(t)/A^-7-)] < 0 (12) 

is the nearest-neighbor (ferromagnetic) coupling strength 
along the Trotter direction. 

In the simulations presented in the main text, the number of 
Trotter slices was fixed to 64. Spin updates were performed 
via Wolff-cluster updates 19^ along the Trotter direction only. 
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(c) Weighted planted, uniform (d) Weighted planted, nonuniform 


FIG. 19. Energy minimization success probabilities of decoding using simulated annealing. 
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FIG. 20. Probability of broken qubits for ME (left) and QAC-ME (right) at different penalty values, for the planted and embeddable planted 
instances. Top row: uniform. Bottom row: nonuniform. The dependence on the clause density is monotonic in the uniform case and exhibits 
a maximum in the nonuniform case. 
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FIG. 21. Probability of broken qubits for ME (left) and QAC-ME (right) at different penalty values, for the embeddable planted instances. Top 
row: uniform. Bottom row: nonuniform. The dependence on the clause density is monotonic in the uniform case and exhibits a maximum in 
the nonuniform case. 
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FIG. 22. Probability of broken qubits for ME (left) and QAC-ME (right) at different penalty values, for the weighted planted and deformed 
embeddable instances. Top row: uniform. Bottom row: nonuniform. The dependence on the clause density is monotonic in the uniform case 
and exhibits a maximum in the nonuniform case. 
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FIG. 23. Probability of broken qubits for ME (left) and QAC-ME (right) at different penalty values, for the deformed embeddable instances. 
Top row: uniform. Bottom row: nonuniform. The dependence on the clause density is monotonic in the uniform case and exhibits a maximum 
in the nonuniform case. 
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PIG. 24. Range and frustration as a function of the clause density. Note the correlation of frustration (range) to the number of broken qubits 
obtained with the nonuniform (uniform) encoding choices of Eigs. |20|2^ 
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FIG. 25. The four panels show the optimal penalty strength 7 for each set of instances and encoding choice (uniform/nonuniform). The optimal 
7 is fairly constant over the different values of the clause density a. 
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(a) Embeddable planted, uniform 


(b) Embeddable planted, nonuniform 




(c) Embeddable planted, uniform. 


(d) Embeddable planted, nonuniform. 


FIG. 26. Comparison between the decoding strategies for the embeddable planted instances. 
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FIG. 27. Comparison between the decoding strategies for the deformed embeddable instances. 



FIG. 28. Mean size of domains connected to a central spin for various N x N x 2 2LG graph sizes for the worst-case percolation probability 
= 0.375. The mean size of connected domains saturates at a size of about 25 spins. We used 10^ random configurations per point. 
































